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Abstract. This paper is part of a series reporting the results of the ATESP radio survey obtained at 1.4 GHz with the Australia 
Telescope Compact Array. The survey consists of 16 radio mosaics with ~ 8" x 14" resolution and uniform sensitivity (lo" 
noise level ~79 /jJy) over the whole area of the ESO Slice Project redshift survey (~ 26 sq. degrees at 8 ~ —40°). The ATESP 
survey has produced a catalogue of 2960 radio sources down to a flux limit (6fJ) of ~ 0.5 mjy. 

In this paper we present the 1 .4 GHz log A' — log S relation derived from the ATESP radio source catalogue. The possible causes 
of incompleteness at the faint end of the source counts are extensively discussed and their effects are quantified and corrected 
for. The ATESP counts are consistent with others reported in the literature, even though some significant discrepancies are 
• present at low fluxes (below a few mJy). We investigate whether such discrepancies may be explained in terms of field-to-field 

Oh anisotropies, considering the fact that all the existing sub-mjy surveys cover small areas of sky (from a fraction of square degree 

to a few square degrees). We stress that the ATESP survey, covering 26 sq. degrees, provides the best determination of source 
5_l ■ counts available today in the flux range 0.7 ^ S\ 4 ghz ~ 2 mjy. 
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1 . Introduction some specific flux range. 

Normalized differential source counts derived from deep 

S3 The statistical study of radio sources allows us to address a L4 GHz suryeys show a flattening below a few mJy xhis 

number of astrophysical and cosmological problems, ranging chang£ fa slop£ h usually interpreted as ±e result of the 

from the evolution of classical radio galaxies and quasars to emergence of a new population which does not appear at 

the properties of starburst galaxies as a function of cosmic higher flux densitieSi where the counts are dominated by the 

time. The radio source counts are the most immediate product dassical powerM radio galaxies and quasars (99% above 60 

which can be derived from a radio survey and reflect the ^ windhorst et al 1990) xhis new population is thought to 

statistical properties of the radio source populations. The be a mixture of several different types of objects (faint AGNs, 

slope of the counts is determined essentially by the relative normal spirals and ellipticaK starburst galaxies , e tc.) whose 

contribution of different types of sources at every flux, which rdative importance changes as a function of me flux itself . 

is the result of their luminosity functions at various redshifts. Unfortunately the existing samples of faint radio sources ^ 

The source counts therefore represent the most immediate confined tQ regions of relatively small angular sizes and no 

observational constraint to evolutionary models of radio complete optical f oUow _ up of such samples has been carried 

sources. An accurate determination of the source counts, their QUt y{Jt (see £ g Gregorini & Prandoni 20 00 for a review). It 

slope and normalization, is however necessary in order to foUows ^ conclusions about the composition of the faint 

make this constraint really useful. The source counts, together radio SQUrces are bafjed Qn smaU and Qnly partially identified 

with measured local radio luminosity functions, can then be sampleSi biased fey ^ fect ^ Qnly ^ brightest optical 

used to make predictions (about the redshift distribution, for counterparts have spectral in f ormation . Many questions are 

instance) which can be verified once an extensive program of therefore gtm open about ±& true nature and evolution of the 

optical identification and spectroscopy has been completed in Mnt radio source population and very little is known about the 
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luminosity properties and redshift distribution of such sources. 
The ATESP 1.4 GHz survey (Prandoni et al. 2000a, paper 
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I) - carried out with the Australia Telescope Compact Array in 
the past years - consists of 16 radio mosaics with ~ 8" x 14" 
resolution and uniform sensitivity (la noise level ranging from 
69 /j Jy to 88 /jJy depending on the radio mosaic) covering a 

narrow strip of 26° x 1° at 5 40°. The ATESP survey has 

produced a catalogue of 3172 radio components, correspond- 
ing to 2960 distinct radio sources, down to a flux limit (6c) 
of - 0.5 mJy (Prandoni et al. 2000b, Paper II). The ATESP 
sample represents the largest sample (26 sq. degr.) sensitive to 
sub-mJy fluxes available so far. It offers a unique opportunity 
to investigate the nature and evolution of the mJy and sub-mJy 
populations and provides a very accurate determination of the 
radio source counts down to sub-mJy fluxes (the most accurate 
in the range 0.7 — 2 mJy). In this respect the ATESP survey 
fully complements the FIRST (White et al. 1997), which 
provides the most accurate determination of the source counts 
between 2 and 30 mJy. 




This paper is organized as follows. The catalogue com- 
pleteness is discussed in Sect. [| The ATESP number counts 
are derived in Sect. || and the possible effect of clustering on 
the faint end of the counts is discussed in Sect. ^. Results are 
summarized in Sect. || 

2. Catalogue Completeness 

The detection of a source depends on the ratio of the measured 
peak flux density]] over the local noise level. In order to use 
the ATESP source catalogue for statistical purposes (like the 
derivation of the source number counts), the possible causes 
of incompleteness in the catalogue have to be recognized and 
taken into account. 

2.1. Visibility Area 

Radio mosaics usually have a varying sensitivity over the re- 
gion surveyed. This means that the area over which a source 
can be detected (visibility area, Katgert et al. 1973) increases 
with the source peak flux. The ATESP survey is characterized 
by a noise distribution which is fairly uniform within each mo- 
saic and from mosaic to mosaic (< 10% variations). Average 
mosaic noise values vary from ~ 69 /jJy to ~ 88 /jJy (see 
Table Bl in paper II) with an average value of ~ 79 /jJy (la). 
Also, the noise can generally be considered Gaussian, with 
only a few exceptions around strong sources (S pea k > 50 — 100 
mJy) where higher noise levels (~ 100 /jJy) can be occasion- 
ally found (see Paper I and II for more details). 
The visibility area of the ATESP survey is shown in Fig. [I]. As 
expected, the fraction of the total area over which an ATESP 
source of a given peak flux density can be detected, increases 
very rapidly between 0.4 and 0.5 mJy and becomes equal to 1 
at Speak > 0.8 mJy. This is a consequence of the ATESP observ- 
ing strategy (uniform sensitivity in the whole region surveyed) 
and of the fact that high noise regions around bright sources do 
not affect the source detection above 0.8 mJy. 

1 Peak flux density is always intended as brightness in units of mJy 
per beam solid angle. 
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Fig. 1. Visibility area of the ATESP survey. Fraction of the total 
area over which a source with given measured peak flux density 
can be detected (i.e. S™^ s > 60). 

2.2. Systematic Effects 

Two additional effects could be responsible of partial incom- 
pleteness at low signal-to-noise ratios: the clean bias and the 
bandwidth smearing. The clean bias has been extensively dis- 
cussed in paper I of this series. It is responsible for both total 
and peak flux density under-estimations of the order of 10- 
20% at the lowest flux levels (5 pea k < 10a) and gradually dis- 
appears going to higher fluxes (no effect for 5 pea k > 50 -100c). 
The bandwidth smearing, extensively discussed in paper II, 
could allow for an extra 5% under-estimation of the peak flux 
densities only. 

In paper II (Appendix B) we gave a formula (Eq. (Bl)) to 
correct the measured (peak and total) fluxes for the clean bias 
and bandwidth smearing. Such a correction varies from source 
to source depending on the radio mosaic and on the source 
signal-to-noise ratio. However, peak flux densities are at most 
underestimated by a factor of ~ 25% due to the combined ef- 
fects of clean bias and smearing (see paper I and II). It follows 
that at most S™ 3 | s ~ 0.75Sp°" k . This means that any source with 
corrected peak flux density J; 6c/0.75 will have effective mea- 
sured peak flux densities still above the detection threshold of 
the ATESP survey. In other words, incompleteness due to clean 
bias and bandwidth smearing does not affect the ATESP cata- 
logue at corrected peak fluxes > 0.55 — 0.70 mJy, the precise 
value depending on the radio mosaic. 

2.3. Resolution Bias 

In deriving the source number counts, the completeness of the 
catalogue in terms of total flux density is very important. Such a 
completeness depends on the source angular size and its com- 
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Fig. 2. Left Panel. Angular size (geometric mean after deconvolution) for the 3 172 ATESP radio components as a function of the 
measured total flux density. The maximum size (0 m ax) that a source of a given S™£jf can have before dropping below the limiting 
flux of the survey is indicated by the dotted line. Also indicated is the minimum angular size (0 m j n ), below which deconvolution 
is not considered significant (solid line). The lines are drawn assuming for the survey a 6o limit of 0.41 mJy, obtained assuming 
the lowest average noise value measured in the ATESP radio mosaics (i.e. o = 69 /jJy, see Table Bl in paper II). All unresolved 
sources are indicated in the figure by dots at © = 0. Right Panel. Same as in the left panel but as a function of the measured peak 
flux density. Here the dotted line indicates the maximum scale (0' ma x) at which the ATESP survey is sensitive due to the lack of 
baselines shorter than 500 m, while the solid line still indicates the minimum angular size (® m i n ), below which deconvolution is 
not considered significant. 



ponent flux ratio. A resolved source of given S to tai will drop 
below the 6o peak flux density cut-off more easily than a point 
source of same Stotai- This is the so-called resolution bias. 
Eq. (1) of Paper II can be used to give an approximate esti- 
mate of the maximum size (@ m ax) a source of given S to tai can 
have before dropping below the 6o limit of the ATESP cata- 
logue. In Fig. H (left panel) we plot the angular size (©) of the 
ATESP radio sources (or source components) as a function of 
the measured total flux density. © is defined as the geometric 
mean of the (Gaussian) major and minor deconvolved axes of 
the ATESP sources as listed in the ATESP catalogue (see Paper 
II). We notice that for unresolved sources (dots at © = 0) we as- 
sume Stotai = Speak- As expected, the angular sizes of the largest 
ATESP sources are in good agreement with the estimated © max 
(dotted line). 

The @ m ax — S to tai constraint is complemented by a second one 
which is a function of the peak flux and takes into account, 
through the so-called visibility function, the maximum scale at 
which the ATESP survey is sensitive due to the lack of base- 
lines shorter than 500 m (see paper I, Sect. 4, for more details). 
This latter constraint (®' ma x) is indicated by the dotted line in 
the right panel of Fig. || and, as expected, follows the angular 
size distribution of the largest ATESP sources as a function of 
their measured peak fluxes. The @' max — S pea k relation is im- 
portant at high fluxes 20 mJy), where it provides a more 



stringent constraint than the ® ma x — S to tai relation. 
A third constraint on source angular sizes should be consid- 
ered when dealing with resolution incompleteness effects. As 
discussed in Paper II, the deconvolution efficiency depends on 
the size of the ATESP survey synthesized beam (~ 8" x 14") 
and increases with the source signal-to-noise ratio: at the low- 
est peak flux levels only sources with © 8" — 10" are reli- 
ably resolved; on the contrary, at the highest peak flux den- 
sities, sources with angular sizes as small as ~ 2" can be 
resolved. Eqs. (1) and (2) of Paper II can be used to derive an 
approximate estimate of the minimum angular size (@ m i n ) that 
is reliably resolved as a function of the source peak flux. © m ; n is 
indicated by the solid lines in Fig. || (since Stotai = S pea k for un- 
resolved sources, we have plotted it in both panels). The © mm 
constraint is important at low flux levels where @ ma x and 0' max 
become both unphysical (i.e. — > 0). 
The three constraints discussed above (0 max , 0'max 

and ©min) 

have been used to define an overall angular size upper limit, 
©lim, as a function of the ATESP source flux density: 

©lim = max[@min,min(@max,© / max)] ■ (1) 

In order to estimate the incompleteness in the ATESP catalogue 
due to resolution bias, we need to make an assumption about 
the true angular size distribution of radio sources as a func- 
tion of flux. Following Windhorst et al. (1990), we have as- 
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Fig. 3. Fraction of sources with angular size larger than 
©lim (h{> ©lim)) as a function of the total source flux 
Stotai (computed as proposed by Windhorst et al. 1990). 
©lim represents the maximum size a source of given 5 tota i 
can have in order to be detected in the ATESP survey. 
Different curves are found for different definitions of ©i; m (see 
text): ©i im = max[@ min ,min(@ max ,0 / m ax)] (solid line); ©i im = 
max[2",min(© max ,©' max )] (dashed line); @i im = max (dotted 
line). Also shown is the curve obtained by assuming ® me d = 2" 
for sub-mJy sources, as discussed in the text (dot-dashed line). 
The inner panel shows the effect of the two different defini- 
tions of me d on the form of the resolution bias correction to 
be applied to the source counts. 



0.30 
1.4 GHz 



sumed h{> @ ]im ) = e-'" 2 (<W®m«i) 62 and © med = 2 " . $ 
(S in mJy). @ me d represents the estimated source median an- 
gular size at a given flux density. The Windhorst et al. equa- 
tions together with Eq. ([[]) allowed us to estimate the fraction 
of sources of given S to ta\ with angular sizes larger than ©i; m and 
thus lost by the ATESP survey (solid line in Fig. The rising 
slope of the curve at high flux densities is due to the constraint 
provided by & max - On the other hand, the decreasing trend at 
low fluxes is a consequence of the looser constraint provided 
by ©min- In fact, introducing © m ; n in the equation takes into ac- 
count the effect of having a finite synthesized beam size (that is 
©lim > > at the survey limit) and a deconvolution efficiency 
which varies with the source peak flux. 

In determining the resolution bias effect, a correct definition 
of @i; m is very important. The dashed and the dotted curves in 
Fig. |5] demonstrate how critical is the determination of the res- 
olution bias effect at low signal-to-noise ratios. Slightly chang- 
ing the definition of @i; m can in fact produce very different 
results. The dashed line shows the effect on h{> ©ii m ) when 
a fixed value for © m ; n in Eq. ([[]) in assumed. In this case 
©min was set equal to 2", that is the minimum angular size 
that is reliably deconvolved at the highest flux densities, in the 



case of the ATESP survey. The dotted line shows, on the con- 
trary, the effect of letting ©i; m — > for S — > 6o, by defining 
©lim = ©max (we also neglect the effect of @' max ). The fraction 
of lost sources will approach 1 at the survey limit (dotted line), 
while it monotonically decreases going to higher fluxes. 
We have also investigated the effect of assuming a constant 
value © me d = 2" for sources with S < 1 mJy, as suggested 
by recent studies (e.g. Windhorst et al. 1993, Richards 2000). 
The corresponding h(> @H m ) is indicated by the dot-dashed 
line in Fig. [| As shown in the figure, this more realistic as- 
sumption results in a larger fraction of lost sources at S < 1 
mJy (10% instead of 6% at the survey limit). The uncertainty 
in the form of © me d at sub-mJy fluxes obviously translates in 
an uncertainty in the form of the resolution bias correction 
(c = 1/[1 — h(> ©iim)]) to be applied to the source counts (see 
inner panel of Fig. |J). At sub-mJy fluxes, this upward system- 
atic uncertainty (< 4%) has to be (quadratically) added to a 
more general 10% uncertainty in the resolution bias correction 
(as indicatively quantified by Windhorst et al. 1990). 

3. The ATESP Source Counts 

We used the 6o ATESP catalogue (Paper II) to derive the 
differential ATESP source counts down to 0.41 mJy. In 
computing the counts we have used the integrated flux for 
extended sources and the peak flux for point-like sources. Each 
source has been weighted for the reciprocal of its visibility 
area (£2(S pea k)/i2 tota i, see Fig. [l]), that is the area over which 
the source could be detected. We notice that for tSp ea k 

> 0.8 

mJy a source can be counted over the whole survey area. 

In Fig. ^ we present the ATESP counts derived from the 
3172 ATESP components (empty circles) and from the 2960 
distinct radio sources defined from them, that is counting 
any multi-component source as one (filled circles). In this 
figure no corrections were applied to take into account either 
the s yste matic under-estimations of th e flu x densities (see 
Sect. ~L2) or the resolution bias (see Sect. |2.3| ). Horizontal bars 
represent la errors in the measured flux (o(5 meas )), calculated 
as discussed in Paper II. Vertical bars represent Poissonian 
errors (calculated following Regener 1951). 
In the ATESP catalogue a source is defined as multiple when- 
ever two or more components are found closer than 45" (even 
though some exceptions are present, see Paper II for more 
details). Applying this distance constraint we expect ~ 20% 
contamination by random superpositions. Fig. ^ shows how 
possible spurious multiple sources could affect the ATESP 
source counts. 

In Fig. |] we show the effect of the flux and resolution bias 
corrections on the ATESP source counts (multiple sources are 
again counted as one). Filled circles and triangles indicate 
the counts respectively before and after correcting the flux 
densities for both the clean bias and the bandwidth smearing 
(fluxes are corrected as proposed in Paper II). In the latter case 
horizontal bars represent a(S con ), obtained by propagating 
the errors of the measured quantities S'" eas , k, a and b which 
appear in Eq. (Bl) of Paper II. We note that the plotted values 
of a(S"" T ) have to be considered as average estimates since 
they have been derived in the case of intermediate clean bias 
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Fig. 4. Normalized 1 .4 GHz differential 
source counts derived from the ATESP 
catalogue. Filled circles refer to the 
ATESP radio source catalogue (2960 
distinct objects). Empty circles refer to 
the list of all the ATESP components 
(3172 entries). Horizontal bars repre- 
sent la errors in the measured total 
flux. Vertical bars represent Poissonian 
errors. 



1000 



100 



• ATESP RS (no corr.) 

a ATESP RS 

(c.b. + sm. corr.) 

■ ATESP RS 

(c.b. + sm. + r.b. corr.) 



10 



0.1 



* 



** 





10 
S (mJy) 



100 



1000 



Fig. 5. Normalized 1 .4 GHz differential 
source counts derived from the ATESP 
catalogue (multiple sources are counted 
as one object). The source counts are 
plotted before (filled circles, same as 
shown in Fig. |j) and after correcting 
for clean bias and bandwidth smearing 
(filled triangles). In the latter case hor- 
izontal bars represent la errors in the 
corrected total flux. Filled squares rep- 
resent the counts obtained once the res- 
olution bias is also taken into account. 
Vertical bars include here the uncer- 
tainty in the resolution bias correction. 
In all cases the counts have been de- 
rived down to 0.70 mJy (instead of 0.41 
mJy), in order to avoid the incomplete- 
ness effects introduced by the flux cor- 
rections. 



correction parameters (a = 0.13 ± 0.03 and b = 0.75 ± 0.06, 
see discussion in Sect. 5.3 of Paper I). We also notice that 
a(k) = 0.006. Filled squares in Fig. || show the effect of adding 
the resolution bias correction (computed following Windhorst 
et al. 1990, as discussed in Sect. E3b. In the estimation of 



the vertical bars we have here included the errors due t o th e 
uncertainty in the resolution bias correction (see Sect. 2.3). 
We notice that the ATESP counts presented in this figure 
are derived down to a limit of 0.70 mJy (instead of 0.41 
mJy) in order to prevent the corrected counts to suffer from 
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Fig. 6. Normalized 1.4 GHz differen- 
tial source counts for different sam- 
ples: the FIRST (triangles, White et 
al. 1997), the Marano Field (squares, 
Gruppioni et al. 1997), the ELAIS 
North and South (stars and crosses, 
Ciliegi et al. 1999 and Gruppioni et al. 
1999), the Phoenix Deep Survey (dots 
and circled dots, Hopkins et al. 1998), 
the Hubble Deep Field + Flanking 
Fields (diamonds, Richards 2000) and 
the Lockman Hole (empty circles, de 
Ruiter et al. 1997). The fit obtained by 
Windhorst et al. (1990) is also shown 
(solid line). The ATESP source counts 
presented here (filled circles), are cor- 
rected for clean bias, bandwidth smear- 
ing and resolution bias (same as filled 
squares presented in Fig. ||). Horizontal 
bars associated to the ATESP counts 
represent la errors in the corrected to- 
tal flux, while vertical bars include the 
uncertainty in the resolution bias cor- 
rection. 



Table 1. ATESP source counts. 
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incompleteness (see discussion in Sect. 



The final ATESP source counts (those corrected for all the 
effects and indicated by filled squares in Fig. ||) are also listed 
in Table [jj where, for each flux interval (AS), the flux geomet- 
ric mean (< S >), the number of sources detected (Ns) and 
the weighted normalized differential counts (nS 2 5 ) are given. 
Also listed are the errors associated to AS and to the normal- 
ized counts (computed as discussed above). In this table fluxes 
are intended as total corrected flux densities. The final ATESP 
source counts have then been compared with the most updated 
previous determinations at 1.4 GHz (see Fig. ^). Also shown is 
the interpolation determined by Windhorst et al. (1990) from 
a collection of 10,575 radio sources belonging to 24 differ- 
ent surveys at 1 .4 GHz, representing the state of the art at that 
time (solid line). As illustrated in Fig. ^| there is consistency 
between the ATESP counts and those obtained by other re- 
cent surveys, with the exception of the Phoenix Deep Survey 
(PDF and PDFS, Hopkins et al. 1998), whose counts at S > 0.7 
mJy are systematically higher than the ATESP counts (and also 
higher than the counts derived from the other surveys presented 
in the figure). 

We notice that according to the ATESP determination the up- 
turn in the source counts should show up at lower fluxes (5^1 
mJy) than indicated by the Windhorst et al. (1990) fit (S <, 5 
mJy). The ATESP counts are in very good agreement with the 
FIRST counts (White et al. 1997, triangles in Fig. ||), which 
are the most accurate available today over the flux range 2-30 
mJy. Our survey, on the other hand, provides the best determi- 
nation of the counts at fainter fluxes (0.7 < S < 2 mJy), where 
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Fig. 7. Normalized 1 .4 GHz differential source counts for the 
16 ATESP mosaics separately (dotted la Poissonian error 
bars). For consistency, these counts are compared only to the 
ones obtained from the surveys covering areas of the same or- 
der of magnitude (a few sq. degrees): the ELAIS North and 
South (filled circles and triangles) and the PDF (circled dots). 
The ATESP counts are not corrected for clean bias, bandwidth 
smearing and resolution bias (Horizontal error bars have not 
been drawn to prevent overcrowding). 

the FIRST becomes incomplete. The ATESP counts can thus 
provide an useful observational constraint on the evolutionary 
models for the mjy and sub-mJy populations. 

Below a few mjy the number counts presented in Fig. |] 
show a large spread. Such a spread can be due to technical 
reasons, like different corrections for resolution bias applied 
by different surveys (e.g. Hubble Deep Field, Phoenix Deep 
Survey, ATESP) or no correction at all (e.g. ELAIS, Marano 
Field). Note, however, that discrepancies are found even within 
the same survey: the deeper PDFS counts, for instance, are 
not consistent with the PDF counts in the overlapping flux 
range. This example suggests that the effect of field-to-field 
anisotropics and/or clustering could be very important, since 
the faintest samples typically cover very small regions of sky 
(<< 1 sq. degree). This latter effect does not apply to the 
ATESP sample which covers 26 sq. degrees. 

4. Clustering Effects on Deep Radio Counts 

The ATESP survey allows us to test whether field-to-field 
anisotropies and/or clustering effects could explain, at least in 
part, the existing spread in source counts at low fluxes (5^1 
mjy). The largest sub-mJy surveys cover areas of the order of 
a few square degrees (typically 3 sq. degr.), while the deepest 
sub-mJy surveys cover only a fraction of square degree (typi- 
cally 0.33 sq. degr.). We have therefore tested these two typi- 



Fig. 8. Normalized 1.4 GHz differential source counts for 0.33 
sq. degree ATESP sub-regions (dotted la Poissonian error 
bars). For consistency these counts are compared only to the 
ones obtained from the surveys covering areas of the same or- 
der of magnitude (a fraction of sq. degree): the PDFS (circled 
dots), the Marano Field (triangles), the HDF+HFF (diamonds) 
and the Lockman Hole (filled circles). The ATESP counts are 
not corrected for clean bias, bandwidth smearing and resolu- 
tion bias (Horizontal error bars have not been drawn to prevent 
overcrowding). 

cal scales for sub-mJy radio surveys. The larger scale has been 
tested by deriving the counts for each of the 16 ATESP mosaics 
separately (mosaic areas range from 1.3 to 2 sq. degrees, see 
paper I). The smaller scale, on the other hand, has been tested 
by splitting the ATESP survey in 40' x 30' (0.33 sq. degrees) 
sub-regions. Note that no correction for clean bias, bandwidth 
smearing and resolution bias has been applied in this analysis 
in order to derive the ATESP counts down to the detection limit 
of the survey (6a) and to have a more meaningful comparison 
with the fainter surveys. 

The results are reported in Figs. ^ and || (dotted error 
bars). For consistency, the ATESP counts obtained on these 
two scales have been compared only to surveys covering areas 
of the same order of magnitudes: the ELAIS and PDF for the 
larger scale (see Fig. g) ; the MF, LH, PDFS and HDF+HFF 
for the smaller scale (see Fig. ^|). In either case the ATESP 
counts are in good agreement with the ones derived from the 
comparison samples, at least down to the fluxes probed by the 
ATESP. This is particularly true if we consider that no cor- 
rections for clean bias, bandwidth smearing or resolution bias 
have been applied to the ATESP counts. It is therefore probable 
that field-to-field anisotropies do not play an important role in 
these samples. The only exception is represented by the PDF 
counts (circled dots in Fig. 0), which remain systematically in 
excess (in particular at fluxes 0.5 & S & 1 mjy). On the other 
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hand, the excess shown by the PDFS counts (circled dots in 
Fig. ||) cannot be directly probed by the ATESP survey. From 
the analysis of Fig. [7] it seems very hard to explain the PDF 
counts in terms of Poissonian fluctuations. We have applied the 
Kolmogorov-Smirnov 2-sample test to the PDF and the ATESP 
counts shown if Fig. ^ (we have considered the counts derived 
from the 16 ATESP mosaics as a unique larger sample) and we 
found that the probability of the two samples being drawn from 
the same parent population is 0.2%. This leads us to conclude 
that either the PDF sample has a significant excess of sources or 
some other unknown systematic effect is present. It is however 
worth to notice that in general systematic discrepancies like 
the ones shown by the counts derived from the Phoenix Deep 
Survey cannot be explained in terms of different assumptions 
for resolution bias, because such an effect does strongly affect 
the c ount s only at low signal-to-noise ratios (see discussion in 
Sect. Oh. 



5. Summary 

We have presented the logN — log S relation derived from the 
ATESP catalogue. Particular emphasis has been given to the 
possible causes of incompleteness at the faint end of the source 
counts. In particular we have discussed the effect of systematic 
under-estimations in the measured flux densities due to clean 
bias and bandwidth smearing and we have quantified the in- 
completeness of the catalogue in terms of total flux density 
(the so-called resolution bias). Once corrected for such effects, 
the ATESP counts are consistent with the counts derived from 
other recent surveys. The only discrepant counts are those ob- 
tained from the Phoenix Deep Survey (Hopkins et al. 1998), 
which are systematically higher. It is worth to notice that the 
ATESP counts point towards lower fluxes (5^1 mJy) for the 
upturn than those (S &5 mJy) indicated by the Windhorst et al. 
(1990) best fit. 

At fluxes below 1 mJy, the 1.4 GHz number counts de- 
rived from different surveys show a large spread. Such a spread 
can be due to technical reasons, like different corrections for 
resolution bias or no correction at all. The effect of field-to- 
field anisotropies and/or clustering could also play an impor- 
tant role, since the faintest samples typically cover very small 
regions of sky (<< 1 sq. degree). This latter effect does not ap- 
ply to the ATESP sample which covers 26 sq. degrees. We have 
investigated such a possibility by splitting the ATESP catalogue 
in smaller sub-samples. In particular we have tested two typi- 
cal scales for sub-mJy radio surveys: a few sq. degrees for the 
largest surveys (1 — 2 sq. degr. in this case) and ~ 0.33 sq. de- 
grees for the deepest surveys. As a result we have found that on 
both scales the spread in the counts can be simply explained in 
terms of Poissonian fluctuations, at least down to the flux densi- 
ties probed by the ATESP. The only exception is represented by 
the Phoenix Deep Survey whose counts remain systematically 
higher than the ATESP counts. Such a discrepancy is statisti- 
cally relevant and draws us to conclude that clustering and/or 
some other unknown systematic effect should play an impor- 
tant role in the case of the PDF sample. 

The ATESP survey provides the best determination of the 
1 .4 GHz source counts in the flux range 0.7-2 mJy, and in this 



respect it complements the FIRST (White et al. 1997) which 
provides the most accurate counts available at higher fluxes (2- 
30 mJy). The ATESP source counts can thus be used to set ob- 
servational constraints on the evolutionary models for the mJy 
and sub-mJy populations. 
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